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Abstract 

The long-standing model-independent annual modulation effect measured 
by DAMA Collaboration is examined in the context of asymmetric mirror 
dark matter, assuming that dark atoms interact with target nuclei in the de¬ 
tector via kinetic mixing between mirror and ordinary photons, both being 
massless. The relevant ranges for the kinetic mixing parameter are obtained 
taking into account various existing uncertainties in nuclear and particle 
physics quantities as well as characteristic density and velocity distributions 
of dark matter in different halo models. 


Keywords: Scintillation detectors, elementary particle processes. Dark Matter 

PACS numbers: 29.40.Me - Scintillation detectors; 95.30.Cq - Elementary par¬ 
ticle processes; 95.35.-|-d - Dark matter (stellar, interstellar, galactic, and cosmo¬ 
logical) . 

1 Introduction 

Annual-modulation effect, as expected from the relative motion of the Earth with 
respect to the relic particles responsible for the Dark Matter (DM) in the galactic 
halo mm, has been measured by DAMA Collaboration using the highly radiopure 
Nal(Tl) detectors of the former DAMA/Nal O 01 [5l [6l [71 El 13 [13 [HI 03 03 03 
[I3030710303[23[2ll[23[23[23[23[23[271[23[23[33[3ll[33 and of the 
second generation DAMA/LIBRA [33[Ml[33[33[37l[33[33[13[lIl[13[13[13[l5] 
apparata. Measurements lasting for 14 annual cycles with an increasing exposure 
which cumulatively is equivalent to 1.33 ton x year confirm the annual modulation 
effect at a confidence level of 9.3a (SB]- No systematic effect or side reaction 
which could mimic the exploited signature, i.e. which would be able to account 
for the whole modulation amplitude and simultaneously satisfy all of the many 
peculiarities of the signature, has been found by the collaboration itself, neither 
it was suggested by anyone else over more than a decade. 
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The annual-modulation effect measured in DAMA experiments is model- 
independent. In other words, the annual modulation of the event rate is an ex¬ 
perimentally established fact, independent on theoretical interpretations of the 
identity of dark matter and specifics of its interactions. It can be related to a 
variety of interaction mechanisms of DM particles with the detector materials (see 
for example Ref. [101 )• In this paper we limit our analysis to the case where the 
signal is induced by atomic-type dark matter candidates from asymmetric mirror 
sector. 

Nowadays the concept that DM may come from a hidden (or shadow) gauge 
sectors which have particle and interaction content similar to that of ordinary par¬ 
ticles becomes increasingly popular. Such a dark sector may consist of elementary 
leptons (analogues of our electron) and baryons (similar to our proton or neutron) 
composed of shadow quarks which are confined by strong gauge interactions like 
in our QCD. These two types of particles can be combined in atoms by electro¬ 
magnetic forces mediated by dark photons. The stability of the dark proton is 
guaranteed by the conservation law of the related baryon number, as the stabil¬ 
ity of our proton is related to the conservation of the ordinary baryon number. 
On the other hand, the cosmological abundance of DM in the Universe can be 
induced by the violation of such baryon number in the early Universe which could 
produce dark baryon asymmetry by mechanisms similar to those considered for 
the primordial baryogenesis in the observable sector. In this respect, such type of 
DM is also known as asymmetric dark matter. 

Historically, the simplest model of such shadow sector, coined as mirror world, 
was introduced long time ago by reasons not much related to dark matter: it was 
assumed that there exists a hidden sector of particles that it is exactly identical 
to the ordinary particle sector, modulo parity transformation: for our particles 
being left-handed, the parity can be interpreted as a discrete mirror symmetry 
which exchanges them with their right-handed twins from parallel mirror sector 
[ig szi ill 09]. Thus, all ordinary particles described by the Standard Model 
SU{3) X SU{2) X U{1), the electron e, proton p, neutron n, photon 7 , neutrinos 
u etc., should have mass-degenerate invisible twins: e', p', n', 7 ', u' etc. which are 
sterile to our strong, weak and and electromagnetic interactions but have instead 
their own gauge interactions, described by the mirror copy of the Standard Model 
SU{3y X SU{2y x U{iy with exactly the same quantum numbers and coupling 
constants. In this case, the microphysics of the mirror matter should be exactly 
the same as that of ordinary matter, at all levels from particle to atomic physics. 

Naively, then one would expect that mirror world would be identical to the 
observable world also in cosmology, which situation could be disfavored by the 
following simple arguments: 

(a) once mirror baryon asymmetry is generated by the same physics as the 
ordinary one, then the cosmological density of the mirror baryons in the Universe 
was expected to be the same as that of the ordinary baryons, = Qb- This 
would not be sufficient for explaining the whole amount of DM. 

(b) the mirror sector with the same abundances of mirror photons and neutri¬ 
nos as the ordinary one would strongly disagree with the Big Bang Nucleosynthe¬ 
sis (BBN) limits on the effective amount of light degrees of freedom. It would be 
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equivalent to effective amount of extra neutrinos AA^gff = 6.15. 

By these reasons, and also because of self-interacting and dissipative nature of 
the mirror matter, it was not considered as a serious candidate for dark matter 
for a long time. However, as it was shown in [50], all problems can be avoided 
assuming that after inflation the two sectors were heated to different temperatures, 
and the temperature of the mirror sector T' remained less than the ordinary one 
T over all stages of the cosmological evolution. The BBN limits are satisfied 
if T' < Tl2 or so, while for T' < r/4 mirror matter can represent the entire 
fraction of DM since in this case early decoupling of the mirror photons renders 
mirror baryons practically indistinguishable from the canonic cold dark matter 
(CDM) in observational tests related to the large scale structure formation and 
CMB anisotropies [501 HUES]. Interestingly, the condition T' < T can also lead 
to mirror baryon asymmetry bigger than the ordinary one, and ~ 5 Db can be 
naturally obtained in certain co-baryogenesis scenarios [53l[5a[55l|56|. Therefore, 
the mirror matter can be a viable candidate for DM, despite its collisional and 
dissipative nature, and the present situation in fundamental physics and cosmology 
gives new perspectives for testing this intriguing hypothesis with a rich predictive 
power (for reviews, see e.g. [5^ [571 l58] l. 

More in general, the mirror DM can also be presented in the form known 
as asymmetric mirror matter, which assumes that mirror parity is spontaneously 
broken and the electroweak symmetry breaking scale v' in the mirror sector is 
much larger than that in our Standard Model, v = 174 GeV [59l |60l [6T1162]. In 
this case, the mirror world becomes a heavier and deformed copy of our world, 
with mirror particle masses scaled in different ways with respect to the masses 
of the ordinary particles. Taking the mirror weak scale e.g. of the order of 10 
TeV, the mirror electron would become two orders of magnitude heavier than our 
electron while the mirror nucleons p' and n' only about 5 times heavier than the 
ordinary nucleons. Then dark matter would exist in the form of mirror hydrogen 
composed of mirror proton and electron, with mass of about 5 GeV which is a 
rather interesting mass range for dark matter particles [631 . In the context of 

baryo-cogenesis mechanisms [SaiMl 1551156], ~ 5Ds can be simply related to 

the mass ratio ~ 5 between mirror and ordinary baryons [561 [551IM] (see also Ref. 
[65]). Owing to the large mass of mirror electron, mirror atoms should be more 
compact and tightly bound with respect to ordinary atoms. Asymmetric mirror 
model can be considered as a natural benchmark for more generic types of atomic 
dark matter with ad hoc chosen parameters. 

In this paper we discuss the annual modulation observed by DAMA in the 
framework of asymmetric mirror matter, in the light of the very interesting in¬ 
teraction portal which is kinetic mixing of two massless states, ordinary 

photon and mirror photon [661IT9] . This mixing mediates the mirror atom scat¬ 
tering off the ordinary target nuclei in the Nal(Tl) detectors at DAMA/LIBRA 
set-up with the Rutherford-like cross sections. The compactness of mirror atoms 
is important here. On the one hand, if these atoms are sufficiently fuzzy, the elec¬ 
tromagnetic scattering between ordinary and mirror nuclei will not be suppressed 
by the mirror atom form-factor. On the other hand, they should be rather com¬ 
pact to render self-interaction of the mirror matter weak enough - otherwise the 
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structure of galactic halo would be strongly affected. Moreover, if the scattering 
cross-section between mirror atoms is such that a/M r-u 10 — 10 cm^/GeV, 

where M is the dark atom mass, some characteristic problems for canonic CDM 
models as the cusp problem or overly large number of small halos within the local 
group can be avoided [671 EH] • 

The paper is organized as follows. In Sect. 2 we give a brief overview of 
asymmetric mirror dark matter models discussing its direct detection via photon- 
mirror photon kinetic mixing. In Sect. 3 details of the analysis are given, while 
in Sect. 4 we discuss the obtained results. 

2 Asymmetric mirror matter and its direct detection 

2.1 Overview of asymmetric mirror matter 

At the basic level, ordinary sector is described by the Standard Model Gsm = 
SU{3) X SU{2) X U{1) containing quarks q = {u,d) and leptons I = of 

three generations and the Higgs doublet H. Then the dark mirror sector must be 
described by the identical gauge group Ggjyj = SU{3y x 51/(2)' x t/(l)' containing 
mirror quarks and leptons, q' = («', d') and I' = {v', e'), and mirror Higgs H'. The 
most general Lagrangian of such Gsm x G'gjyj theory has the form 

/Itot = C + C' + Tmix (1) 

where the Lagrangians C = Tgauge + /Ivuk + ^Higgs and C! = ^ ^ffiggs 

describing all interactions in the ordinary and mirror sectors can be rendered iden¬ 
tical by imposing mirror parity, a discrete symmetry that exchanges C ^ C . The 
Tmix describes the possible interactions between ordinary and mirror particles as 
e.g. photon-mirror photon kinetic mixing which shall be discussed later. One can 
make a further step extending towards the concepts of supersymmetry (SUSY) 
and grand unification theory (GUT). In view of the gauge coupling unification, 
envisaging that both sectors, ordinary and dark, in the ultraviolet limit are de¬ 
scribed by a SUSY GUT G x G' with two identical gauge factors that can be 
SU{5), 50(10) or 5t/(6). Then both gauge factors G and G' are spontaneously 
broken down to their standard subgroups at the grand unification scale around 
10^® GeV. Below this scale ordinary and mirror sectors are represented by re¬ 
spective Standard Models (or SUSY Standard Models), SU(3) x 51/(2) x U{1) 
and 51/(3)' x SU{2)' x //(!)', with identical particle contents and identical pat¬ 
terns of interaction constants (gauge and Yukawa) 0 In the case when the mirror 
parity remains unbroken at all levels, so called exact mirror model, we have no 
new parameter: ordinary and mirror particles are degenerate in mass, and the 
microphysics of mirror and ordinary sectors is exactly the same. 

However, one can consider another possibility when mirror parity is sponta¬ 
neously broken at the level of electroweak interactions so that the vacuum expec- 

^In the following, we do not consider supersymmetric contributions to dark matter. The SUSY 
DM in the form of neutralinos can be destabilized by introducing a tiny violation of i?-parity, or 
perhaps not so tiny but in specific forms which would not affect the matter stability itself, as e.g. 
in Refs. [6^[70| . 
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tation value (VEV) of mirror Higgs {H') = v' is much larger than that of the 
ordinary Higgs {H) = v ^ 174 GeV [59l [601 |6T] H Such a concept of asymmetric 
mirror matter introduces a new parameter, the ratio C = v'/v between the VEVs, 
while the gauge and Yukawa constants in two sectors essentially remain identical. 
When ^ 1, the mirror sector becomes a deformed heavier copy of the ordinary 
sector, with mirror particle masses scaled in different ways as functions of C with 
respect to the masses of their ordinary partners. Namely, the masses of shadow 
quarks and charged leptons are induced by the Yukawa couplings with the Higgs 
H', and so they scale roughly by a factor with respect of the masses of their 
ordinary twins (e.g. m'^ ud — for the electron, up-quark and down-quark, 

if one neglects the renormalization group running factors for quarks from the scale 
{H') down to (H)), the mirror neutrino masses which emerge from the dimension- 
5 operators H'H' quadratic in Higgs field, scale by a factor with respect 
to ordinary neutrino masses generated by similar operators In addition, 

the ‘Eermi’ constant of mirror weak interactions becomes smaller by a factor 
with respect to our Eermi constant Gp, while the infrared scale of mirror QCD 
changes roughly as A' ~ A with respect to our QCD scale due to the threshold 
effects of heavier mirror quarks. If C < 10^j for the masses of light mirror quarks 
u' and d! we have m'^ ^ < A', and so for the mass of mirror pions we get roughly 
~ Hence, taking e.g. C = 100, shadow electrons become two orders 
of magnitude heavier than our electron, shadow neutrinos 10^ times heavier than 
our neutrinos, but shadow proton and neutron only about 5 times heavier than 
ordinary nucleons, e.g. m'p ~ due to larger by larger QCD 

scale and larger masses of light mirror quarks. Hence, mirror nucleon masses be¬ 
come about 5 GeV which is an interesting mass range for dark matter particles. 
In addition, due to the large mass gap between the masses of mirror neutron and 
proton, m'^ — m'p S> m'g, the neutron becomes unstable even if it is bounded in 
nuclei, and hence the only stable atom in such shadow sector can be the mirror 
hydrogen composed of mirror proton and electron issiEniiii]. Needless to say, 
that mirror hydrogen can be presented in molecular form, while some its fraction 
can exist in ionized state: if dark electron and proton are sufficiently heavy, the 
matter may remain essentially ionized owing to very small recombination cross 
section. 

The baryogenesis in the two sectors, ordinary and mirror, emerges by the same 
mechanism, since the particle physics responsible for baryogenesis is the same in 
the two sectors (coupling constants, CP-violating phases, etc.). The spectrum of 
mirror particle masses is irrelevant for the baryon asymmetry in the mirror sector 
as far as the effective scale of interactions relevant for baryogenesis is much higher 

^Mechanisms of discrete symmetry breaking between ordinary and mirror sectors were dis¬ 
cussed e.g. in Refs. |601161|. Seemingly, there emerges a hierarchy problem between the two 
VEVs v' and v, which however can be turned into positive aspect ameliorating the hierarchy 
problem for the ordinary Higgs, in the context of supersymmetric little Higgs or twin Higgs mod¬ 
els | 58l 1711172| . In these models the Higgs in ordinary sector emerges as a pseudo-Goldstone 
boson associated with the breaking of the accidental global SU (4) symmetry between two Higgs 
systems, and thus the mass and VEV of ordinary Higgs 77 can be smaller, within two orders of 
magnitude or so, than the VEV of the mirror Higgs 77', v' ~ 10 TeV. 
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than the mirror electroweak scale. However, the cosmological conditions at the 
baryogenesis epoch can be different (recall that shadow sector must be colder than 
ordinary one). One can consider two cases: 

1) Separate baryogenesis, when the baryon asymmetry in each sector is gener¬ 

ated independently but by the same mechanism. In this case, in the most naive 
picture when out-of-equilibrium conditions are well satisfied in both sectors, one 
predicts rj = ns/n^ and rj' = must be equal, while ~ <C 1, where 

X = T'/T is the temperature ratio between mirror and ordinary worlds in the early 
Universe. In this case, we have Q'^/Qb — /mj<!)x^ where m'^ is the mass of 

shadow nucleon vs. the mass of ordinary nucleon m^v — 1 GeV. Therefore, if e.g. 
X = 0.5, for Q'^/Qb — 5 we need m)y = 40 GeV. In the context of asymmetric 
rescaling of particle masses, this would occur when ~ 2000 or so, in which case 
for mirror electron mass one would have rUg ~ 1 GeV or so. In this case, we are 
obviously in rather heavy range of dark matter. However, one should remark that 
due to different out-of equilibrium conditions in two sectors situation with rj' ^ r] 
can be also obtained in some specific parameter space [50] . 

2) Co-genesis of baryon and mirror baryon asymmetries via B — L and GP- 
violating processes between the ordinary and mirror particles, e.g. by the terms 

ll'HH' in Tmix which also induce mixing between ordinary (active) and mirror 
(sterile) neutrinos, and which can be mediated by heavy “right-handed” neutrinos 
coupled to both sectors as e.g. [S31IM|- This leptogenesis mechanism predicts 
^'b = nB and so for Q'^/Qb — 5 we need m'j^/mi\f ~ 5, which singles out the 
mass of dark atom of about 5 GeV. (Somewhat different leptogenesis mechanism 
also based on the assumption > niN was suggested in Ref. [65].) This would 
occur when ^ ~ 100 in which case the mirror electron mass is m'^ ~ 50 MeV 

[MllSHKMlEl]. 

Let us remark, however, that in the picture of asymmetric mirror world the 
mirror (sterile) neutrinos could also be a natural candidate of dark matter, namely 
warm dark matter. In fact, assuming a minimal mass normal hierarchy for ordi¬ 
nary neutrino masses, one can approximate to current cosmological data with a sin¬ 
gle mass eigenstate with mass my = 0.06 eV. We have then ~ ^ eV ~ 

6 X 10“^, well below the cosmological limits on the neutrino masses. Since the 
mirror neutrino masses scale roughly as a = 1,2,3, then for 

( ~ 10^ or more we would have ~ (C/100)^ ^ 0-6 keV and 

~ C‘^XyY^my/93 eV ps 0.1(C/100)^(x/ 0.25)^, where Xy = T'yjTy is the tem¬ 
perature ratio between the mirror and ordinary neutrinos in the early Universe 
(which can be somewhat smaller than the photon temperature ratio x = T'/T if 
is large) [601161] . Therefore, for C, — 100 and x ~ 1/4 practically the whole budget 
of dark matter will be saturated by mirror neutrinos with mass ~ 1 keV leaving 
practically no space for mirror baryons. Therefore, in the following we assume 
that the mirror sector is cold enough to leave a space for mirror baryons as well. 
In any case, in what follows, we do not require that mirror baryons provide entire 
amount of DM, but we assume that it provides some fraction / of DM which we 
shall keep as an arbitrary parameter. 

How large this fraction can be depends on the mass spectrum of the mirror 
particle and it is limited by the degree of self-interaction of mirror atoms. The self- 
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scattering cross section of dark matter particles should satisfy the (most conserva¬ 
tive) upper limit a/M < 10“^^ cm^/GeV [671168j . For this, mirror atoms should 
be enough compact, with small enough mirror Bohr radius a' = — cl/C 

where we take m'^ <C m'p and a = 5.3 x 10“® cm is ordinary Bohr radius. In the 
case of mirror hydrogen this condition cannot be easily satisfied. The hydrogen 
atom has a long tail potential which in fact is responsible for molecular hydrogen 
as a bound state of two atoms. Therefore, the scattering cross section is resonantly 
amplified when the energy of scattering E gets close to the energy levels of the 
hydrogen molecule. One can parametrize this cross section as (t{E) = F{E)a'‘^. 
At energies exceeding the atom binding energy, E > Eq = a'^mg/2, one has 
E{E) < 1 , but for E Eq the cross-section strongly increases, and we have 
EiE) 100 or so m- Let us consider e.g. the case m'p = 6 GeV for the 
dark proton mass which can be obtained when C ~ 200. Then at typical virial 
velocities v ~ 300 km/s characteristic for large galactic halos, we have FI ~ 3 
keV, which is comparable to E'q when m'^ = Qme — 100 MeV. However, in this 
case we get a/M ~ a'‘^/m'p ~ 10“^^ cm^/GeV which exceeds by an order of 
magnitude the above upper limit {a}/m^ ~ 3 x 10“^^ cm^/GeV for ordinary hy¬ 
drogen). If we take instead / = 2000, we get m!^ ~ 1 GeV and m!p ~ 20 GeV 
or so. In this case we have FI ~ 10 keV against E'q ~ 30 keV, so E{E) ~ 10 
and a/M ~ lOa'^/m), ~ 3 X 10 cm^/GeV which can be acceptable for large 
halos, with v ~ 300 km/s. Moreover, such a self interaction would smooth out the 
cusp providing a central density of halos of about 10“^ Mq/pc^ ~ 0.3 GeV/cm^ 
[671168j . However, for small halos, with e.g. u ~ 30 km/s, we would have E Eq 
and so E{E) > 100 [73], and thus we would get a/M > 3 x 10“^^ cm^/GeV, again 
larger than the above conservative upper limit. Therefore, in the cases discussed 
above mirror hydrogen cannot provide the entire amount of dark matter. In the 
case C ~ 10^ it can be only a subdominant fraction (e.g. / ~ 0.1) while the rest 
of dark matter can be represented e.g. by sterile mirror neutrinos with masses 
of few keV. In the case (/ ~ 10^ mirror hydrogen could be a dominant fraction, 
say / ~ 0.8 or so. However, it would be not effective for small dark matter halos 
(dwarf galaxies) and thus the latter must be formed essentially by mirror neu¬ 
trinos constituting a small fraction, say 20 per cent or so of the whole amount 
of DM. This transition of DM content from large to small halos may provide an 
interesting explanation to the reduced amount of small substructures and “Too 
big to fail” problems. 

More generically, the low energy theory in both sectors could be represented 
by the Standard Model containing two Higgs doublets Hu and Hd, as e.g. in 
the context of supersymmetry or in the model providing an axion solution to the 
strong GP-problem m- In this case the ratio of two VEVs tan/3' = (H'u)/{H'^) 
in mirror sector could be different from tan/3 = {Hu)/{Hd). In other words, 
up and down quark masses in the mirror sector will scale by different factors 
Cu = {H'u)/{Hu) and Cd = {Hd)/{Hd): we have = Cd"ie,M,-r for charged 

leptons and m'^ sb ~ Cd'cnd,s,b for down-type quarks, but = Cu'mu,c,t for 

up-type quarks and for the neutrino Majorana masses. 

For the time being, tan /3 is a free parameter for the MSSM or for the two 
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Higgs model which formally can range from 1 to abont 100 in ordinary sector. 
The case when tan/3' = tan/3, i.e. Cu = Cd = C; was discussed above. However, if 
tan /3' 7 ^ tan /3, other interesting options can appear. For example, taking ~ 10^ 
and Cd 10^, which correspond to the case e.g. tan/3 ~ 10 and tan/3' ~ 1, we 
would have ~ Cd^e — 500 MeV. But now mirror upper quarks are rather 
light, mC ~ lOOm^ ~ 250 MeV, while mirror down quarks become very heavy 
m'^ ~ lO^rurf ~ 5 GeV. In this case down quarks are unstable against /3-decay even 
when they are bounded in hadrons, and thus the lightest baryon in mirror sector 
should exist in the form of A'"*"*’ = u'u'u' bound state with spin 3/2 and mirror 
electric charge -|-2, rather than in the form of the mirror proton. The mass of such 
baryon will be rescaled roughly as A'/A ~ 4 with respect to ordinary A-resonances, 
M/\ ~ 1.2 GeV. On the other hand, because of charge +2, it must form a helium¬ 
like atom with two electrons, having the mass ~ 6 GeV or so. In this case 
the chemical neutrality of the helium gives an interesting indication for the self¬ 
scattering cross section. Taking into account that the helium atom radius is about 
3/5 of the Bohr radius a' = a/Cd, we obtain oa’A'I^'a — 2 x 10“^^ cm^/GeV in 
excellent agreement with the self-scattering limits. Therefore, such exotic helium 
atom could constitute the entire dark matter in the galaxies and in the whole 
universe. 

There can emerge other interesting possibilities. E.g. for a proper choice of 
tan/3' ~ 2tan/3, mirror proton and mirror neutron can be enough degenerate in 
mass, both could be stable and also form mirror nuclei also with rather large 
atomic numbers. Alternatively, if tan /3' > 2 tan /3 or so, mirror proton could 
become heavier than mirror neutron and it can decay into the latter. In this case 
mirror world will be represented solely by mirror neutrons with mass say 5 GeV 
and with strong self-interaction cross-section, an'n' ~ 10“^^ cm^ or so. Hence, we 
have rate cr/M ~ 10“^^ cm^/GeV, which is in perfect agreement with all limits 
on self-scattering DM and yet can solve the problems of galactic halo shape. 

2.2 Photon-mirror photon kinetic mixing portal 

The mirror matter can interact with ordinary matter through the photon-mirror 
photon kinetic mixing. In the context of Gsm x kinetic mixing between 

gauge bosons of two abelian factors, | is allowed. After the electroweak 

symmetry breaking, it transforms into photon ~ mirror photon kinetic mixing term 

(2) 

with e = e cos 9w cos . This mixing does not induce oscillation between ordinary 
and mirror photons as far as both are massless. However, it in fact makes mirror 
particles mini-charged with respect to ordinary electromagnetic interactions: the 
mirror particles acquire electric charges eq. Generically, in this case dimensionless 
parameter e could be order 1. 

On the other hand, there are stringent experimental constraints on the param¬ 
eter e which depend on the masses of mini-charged particles [751 Eg [771 [Tg [791 [80]. 
As we see below, our results for dark matter detection are compatible with the 
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existing limits on the dark matter particle mini-charges. In particular, regarding 
the mass range of mirror particles we are interested in, taking mirror electron 
mass e.g. ~ 100 MeV, we have a direct experimental limit e < 6 x 10“^ 
while cosmological BBN limit from e'’'e s-channel annihilation into mirror 
electron positron pair gives e < 10“^(mg/l GeV)^/^ or so [80] , 

There emerges a question (or naturalness problem), i.e. why e is so small, 
which can be easily resolved by considering a grand unified version of mirror mat¬ 
ter, G X G' with e.g. G = SU{5) or SU{6). In this case, no kinetic mixing is 
possible between the non-abelian gauge bosons of two sectors without the GUT 
symmetry breaking. Hence, the term ([2|) can emerge only from the higher dimen¬ 
sional operator {TjG^'^){Ti'G'^^) where G^'^ are gauge fields of the GUT, S, S' 
are Higgs fields in adjoint representations which break the GUT symmetry down 
to the Standard Model, and Mp is some natural cutoff scale, of order of Planck 
mass or so. Therefore, photo-mirror photon kinetic mixing can emerge only after 
breaking of the GUT symmetries by the VEVs (S) = (S') ~ GeV, and so 

one can naturally have e < 10“^ even if constant k is order of one |56l I58j . 

The basic mechanism in the mirror atom scattering off the ordinary target is 
the following: the mirror nuclei A/"', with a mirror electric charge Z', interacts 
with M ordinary nuclei with the electric charge Z: 

+ + M (3) 

via kinetic mixing of two photons as shown in Fig. [T]with Rutherford-like cross 
section suppressed by the free parameter e <C 1. (The direct scattering of the 
mirror nuclei scattering off the electrons will be irrelevant by obvious reasons). 
The mirror nuclei can be simply a mirror proton, Z' = 1, in the case when mirror 
world is dominated by hydrogen-like atoms. In other situation discussed above, 
when mirror matter is presented by A-atoms, N' can be A'’*”*' particle with Z' = 2. 

The effect of the e' screening will be negligible if the mirror atom is not too 
compact, i.e. the inverse radius of the mirror atom 1/a' ~ am'^, is smaller than the 
transfer momentum q = where Ma is the mass of target atom and Fp is 

recoil energy. In particular, for Na target in DAMA, considering that the relevant 
recoil energy range is 2-6 keV electron equivalent (keVee) which corresponds to 
Efi ~ 6 —20 keV when one takes into account the quenching factor, we have q > 20 
MeV. Thus, for m'^ < 1 GeV the condition 1/q < a' is safely satisfied, and as a 
consequence, the cross section depends on the recoil energy, Up, and the relative 
velocity, v, as da/dEji ~ l/(Fpu)^. 

Gertainly, one could consider also contact interactions between ordinary and 
mirror nucleons due to interaction terms like -^q'y^qq''^^q' between ordinary and 
mirror quarks which can be mediated by extra gauge bosons connecting two sec¬ 
tors, as e.g. of flavor symmetry [83]. However, if these interactions have cross- 
sections large enough for the dark matter direct detection, then the same inter¬ 
actions would bring two sectors into equilibrium in the early universe, violating 
the BBN limit and even overdosing the universe due to large masses of the sterile 

®In the case of exact mirror model, this mechanism was discussed by Foot [smu]. 
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Figure 1: Mirror Nucleus-Nucleus interaction through Photon-Mirror photon por¬ 
tal. 


mirror neutrinos. This can be avoided only by assuming a very low reheating tem¬ 
perature after the inflation. As for the photon-mirror photon kinetic mixing portal, 
the cosmological constraints are easily satisfied without any exotic assumptions 
on the post-inflationary reheating of the universe. By the similar reasons, this 
interaction portal avoids the experimental limits on the dark matter production 
at the LHC and other accelerators. 


2.3 Interaction Rate 


The low-energy differential cross-section of the interaction between mirror and 
ordinary nuclei has the Rutherford-like form: 

dcrA,A' _ Ca,A' ..X 

dEn ~ Elv^ ^ ’ 

where E^ is the energy of the ordinary nucleus recoil, v = |r7| is the relative 
velocity between the mirror nucleus and the ordinary one, and; 


Ca,A' 


2TTe‘^a‘^Z^Z'‘^ 

'^A 


772 772 
J- J^J- 


(5) 


where a is the fine structure constant, Z and Z' are the charge numbers of the 
ordinary and mirror nuclei, Ma is the mass of the ordinary nucleus, and ExifEx) 
(X = A, A') are the Form-factors of ordinary and mirror nuclei, which depend on 
the momentum transfer, q, and on the radius of X nucleus. 

As a consequence the differential interaction rate of mirror nuclei on a target 
composed by more than one kind of nucleus is: 


dU 

dEn 


Naua' 

A,A' 


da A, A' 
dEji 


EA'iv,VE)vd^v = 


( 6 ) 


10 



J.J- Ca,a' f ^A'{v,ve) 

= 2^^AnA'-^ - 

^ J\sJ FI (Er) ^ 

where; i) Na is the number of the target atoms of specie A per kg of detector; ii) 
n-A' = PdmfA'/^A' with pdm local dark matter density, fA' fraction of the specie 
A' in the dark halo, and Ma' mirror nucleus mass; iii) the sum is performed over 
the mirror nuclei involved in the interactions {A') and over the target nuclei in the 
detector {A)] iv) Fa/{v,ve) is the velocity distribution of the A' mirror nuclei in 
the laboratory frame, which depends on the velocity of the Earth in the galactic 
frame; ve- The lower velocity limit Vmin{EE) is 


e) 


\Ma + Ma'YEr 


2MaMI 


A' 


(7) 


As mentioned above, in our model we consider just one specie of mirror nuclei. Our 
benchmark model is the mirror hydrogen [A' = Z' = 1), with mass Ma' — 5 GeV 
and fraction fA' = /• Alternatively, one can consider the helium like A-atom that 
we have discussed above, with A' = \, Z' = 2 and with mass again Ma' — 5 GeV. 
All numerical results, presented below in the case of mirror hydrogen in terms of 
would be equivalent to Z'ef^^"^ in the case of mirror nuclei with Z' > \ with 
the same mass. So, for A-atom one just puts Z' = 2. 

In order to compare the theoretical differential rate with the experimental data, 
one has to take into account the detector response by means of the 1Ca{E\Ee) 
kernel for each A nucleus in the target-detector; E is the detected energy in keV 
electron equivalent (generally in literature indicated simply as keV). Thus, the 
theoretical differential counting rate can be written as; 

g = E/ 'CAElEA§^dE,. (8) 


where is the differential interaction rate on the A nucleus in the detector. The 
ICa{E\Ee) kernel accounts for the detector’s energy resolution (generally through 
a gaussian convolution) and for the transformation of the nuclear recoil energy in 
keV electron equivalent through the use of a quenching factor. For a discussion 
about the quenching factors see later. In particular, that kernel can be written as; 


ICa{E\Er) = / g{E\E')QAiE'\ER)dE' 


(9) 


where; 

g{E\E') = 

takes into account the energy resolution a (generally function of E') by a gaussian 
behaviour, and Qa{E'\Er) takes into account the energy transformation through 
the quenching factor (see later). For example, the latter kernel can be written in 
the simplest case of a constant quenching factor qA as; Qa{E'\Er) = 5{E'—qaEr). 

The expected differential rate depends on the Earth’s velocity in the galactic 
frame, ve, which depends on the time of the year. Projecting VE{t) on the galactic 
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plane, one can write: VE{t) = VQ+v^cos'ycosu}{t — to); here Vq is the Sun’s velocity 
with respect to the galactic halo (n© ~ uq + 12 km/s and vq is the local velocity), 
n© ~ 30 km/s is the Earth’s orbital velocity around the Sun on a plane with 
inclination 7 = 60" with respect to the galactic plane. Furthermore, uj= 27r/T 
with T=1 year and roughly to ~ June 2"''’* (when the Earth’s speed in the galactic 
halo is at maximum). The Earth’s velocity can be conveniently expressed in unit 
of Vo- rj{t) = VE{t)/vo = 7 o + Ar]cosuj{t — to), where - depending on the assumed 
value of the local velocity ranging between 170 and 270 km/s - ?7o=l-04-1.07 is the 
yearly average of rj and Arj = 0.05-0.09. Since Ary <C tjo, the expected counting 
rate can be expressed by the first order Taylor expansion; 


dTZ, , ,, dTZ, , d (d'R\ 


Arycosa;(t — to). 


v=vo 


Averaging this expression in a given energy interval one obtains: 


( 10 ) 


S[r]{t)] = 5[ryo] + 



Ar]COSLv(t — to) = So + SmCOSUj{t — to), 


( 11 ) 


with the contribution from the highest order terms less than 0 . 1 %; Sm and So 
are the modulated and the unmodulated part of the expected differential counting 
rate, respectively. 


3 Details of the analysis 

The data analysis in the Mirror DM model framework considered here allows the 
determination of the ^/fe parameter. As mentioned this corollary analysis is model 
dependent; thus, it is important to point out at least the main topics which enter 
in the y/fe determination and the related uncertainties. In the following the main 
ones are addressed. 

3.1 Phase-space distribution functions of the dark halo 

In order to derive the y/fe parameter a specific phase-space Distribution Func¬ 
tion (DF) of the Mirror Dark Matter in the galactic halo has to be adopted. A 
large number of possibilities is available in literature; thus, this introduces large 
uncertainties in the predicted theoretical rate. In addition, it is also possible the 
presence of non-virialized components, as streams in the dark halo coming from 
external sources with respect to our galaxy [Mj; these latter possibilities are not 
included in the present analysis. In conclusion, it is strongly limiting/arbitrary 
to just consider an isothermal profil^ with local parameters vq = 220 km/s and 
— 0.3GeV/cm^ without taking in consideration other existing possibilities 

“^It is also worth noting that the isothermal halo is an unphysical model; for example, the 
mass would diverge and one has to adopt a by-hand cut-off. Let us remark, however, that fiat 
density profile for the Galaxy within the radius of 10 kpc can be obtained if the DM particles 
have self-interaction cross-section a/M ~ — 10“^® cm^/GeV |671168| . 
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in the distribution of velocity and spatial coordinates permitted by astrophysical 
observations. 

In this paper, we will consider a large (but not exhaustive) class of dark halo 
models, as already done in previous analyses for other DM candidates miiiiiiisi 
ng ttzi [HI da [201 EH Ezi HD] ; these models are summarized in Tabled An extensive 
discussion about some of the more credited realistic halo models has been reported 
in Ref. |lllll4j . In particular, the considered classes of halo models correspond to: 
i) spherically symmetric matter density with isotropic velocity dispersion (A); ii) 
spherically symmetric matter density with non-isotropic velocity dispersion (B); 
iii) axisymmetric models (C); iv) triaxial models (D); v) moreover, in the case of 
axisymmetric models it is possible to include either an halo co-rotation or an halo 
counter-rotation. 

In our analysis we also consider the physical ranges of the main halo param¬ 
eters: the local total DM density pdm and the local velocity vq as previously 
discussed in Ref. m- In particular the range of the possible vq value is from 
170km/s to 270km/s. We will consider for pdmj hs minimal value Pd^ 
maximal one compatible with the given vq. The PTnf 

fined (as in Ref. ffH) as the values, associated to a specific pdm, which provide 
a visible mass maximal and minimal contribution, respectively, to the total mass 
of the dark halo compatible with the astrophysical observations and constraints. 
The particular values for Pd^^ related to the particular DF and the par¬ 

ticular vq considered. See Table III of Ref. m for maximal densities at given 
Vq and given model; for vo= 170 km/s pdm ranges from 0.17 to 0.67 GeV cm“^, 
while for vo= 220 km/s pdm ranges from 0.29 to 1.11 GeV cm“^, and for vo= 
270 km/s Pdm ranges from 0.45 to 1.68 GeV cm“^, depending on the halo model. 
Finally, we consider a DM escape velocity from the galactic gravitational potential 
Vesc = 650 km/s, as often considered in literature; however, it is also affected by 
significant uncertainty. No sizeable differences are observed in the final results 
when a different value of Vesc = 550 km/s is considered. 

3.2 Nuclei and Dark Matter Form factors 

Other important items for the determination of the expected signal counting rate 
are the nuclei and DM form factors. Usually a Helm form factor [851 [56] is con¬ 
sider ecd for each X ordinary and mirror nucleus: 

Hxiqrx) = ^^ 2 ) 

qrx 

where q = (2MxFij)^/^ is the momentum transfer, rx is the effective nuclear 
radius (the normalization in natural units h = c = 1 is understood), s is the nu¬ 
clear surface thickness, and ji is the Bessel function of order 1. This analytical 
expression is sufficiently good for our purposes, especially comparing the uncer¬ 
tainties coming e.g. from the astrophysical side. We consider here s ~ 1 fm, 

®It should be noted that the Helm form factor is the less favorable one e.g. for iodine and re¬ 
quires larger SI cross-sections for a given signal rate; in case other form factor profiles, considered 
in the literature, would be used |14| . the allowed parameters space would extend. 
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Table 1: Summary of the considered consistent halo models [111 I14j . The labels 
in the first column identify the models. In the third column the values of the 
related considered parameters are reported mm-, other choices are also possible 
as well as other halo models. The models of the Class C have also been considered 
including possible co-rotation and counter-rotation of the dark halo. 


Class A: spherical pdm? isotropic velocity dispersion 

AO 

Isothermal Sphere 


A1 

Evans’ logarithmic 

Rc = 5 kpc 

A2 

Evans’ power-law 

Rc = 16 kpc, j3 = 0.7 

A3 

Evans’ power-law 

Rc = 2 kpc, (3 = —0.1 

A4 

Jaffe 

a = 1, /3 = 4, 7 = 2, o = 160 kpc 

A5 

NEW 

a = 1, /3 = 3, 7 = 1, a = 20 kpc 

A6 

Moore et al. 

a = 1.5, /3 = 3, 7 = 1.5, a = 28 kpc 

A7 

Kravtsov et al. 

a = 2, /3 = 3, 7 = 0.4, a = 10 kpc 

Class B: spherical pdm? non—isotropic velocity dispersion 

(Osipkov—Merrit, /3o = 0.4) 


B1 

Evans’ logarithmic 

Rc = 5 kpc 

B2 

Evans’ power-law 

Rc = 16 kpc, /3 = 0.7 

B3 

Evans’ power-law 

Rc = 2 kpc, f3 = —0.1 

B4 

Jaffe 

a = 1, /3 = 4, 7 = 2, a = 160 kpc 

B5 

NEW 

a = 1, /? = 3, 7 = 1, a = 20 kpc 

B6 

Moore et al. 

a = 1.5, P = 3, j = 1.5, a = 28 kpc 

B7 

Kravtsov et al. 

a = 2, /3 = 3, 7 = 0.4, a = 10 kpc 

Class C: Axisymmetric pdm 

Cl 

Evans’ logarithmic 

1—1 

II 

CD 

II 

C2 

Evans’ logarithmic 

Rc = 5 kpc, q = I/V2 

C3 

Evans’ power-law 

Rc = 16 kpc, q = 0.95, /? = 0.9 

C4 

Evans’ power-law 

Rc = 2 kpc, q = I/V 2 , /3 = -0.1 

Class D: Triaxial pdm (q = 0.8, p = 0.9) 

D1 

Earth on maj. axis, rad. anis. 

6 = -1.78 

D2 

Earth on maj. axis, tang. anis. 

J = 16 

D3 

Earth on interm, axis, rad. anis. 

J = -1.78 

D4 

Earth on interm, axis, tang. anis. 

(5 = 16 


rx = Y^rg — and rg = 1.2^4^/^ fm; in case of light nuclei, as e.g. the mirror 
Hydrogen is, we use s ~ 0.9 fm and rx = fm. However, let us note that 

the mirror proton is even more compact that ordinary proton and thus we can 
safely take the nuclear form factor equal to one. In the analysis some uncertainties 
on the nuclear radius and on the nuclear surface thickness parameters in the Helm 
SI form factors have been included (see e.g. naEZj). 
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3.3 Quenching factors and Channeling effect 

In the present analysis, we consider with the appropriate care the uncertainties in 
the quenching factors. A precise experimental determinations of these quantities 
are difficult for all kind of detectors. In fact, generally the direct measurements 
of quenching factors are performed with reference detectors, and - in some cases 
- with reference detectors having features quite different from the ones used in 
the underground running conditions; in other cases the quenching factors are not 
even measured at all. In addition it should be noted that the quenching factor 
value is a feature of each specific detector and not a fixed property of a given 
material. Moreover, the real nature of these measurements and the used neutron 
beam/sources may not point out all the possible contributions or instead may cause 
uncertainties because e.g. of the presence of spurious effects due to interactions 
with dead materials as e.g. housing or cryogenic assembling, if any; therefore, they 
are intrinsically more uncertain than generally derived. A discussion dedicated to 
the case of Na and I quenching factors in DAMA experiments has been given in 
section II of Ref. [37] ; analogous / similar discussions should be pursued for every 
other case. In fact, the related uncertainties affect all the results both in terms of 
exclusion plots and in terms of allowed regions/volumes; thus, comparisons with a 
fixed set of assumptions and parameters values are intrinsically strongly uncertain. 

According to Ref. [37112II, in the present analysis three possibilities for the Na 
and I quenching factors have been considered: Q/) the quenching factors of Na 
and I “constants” with respect to the recoil energy Er: q^a — 0.3 and qi ~ 0.09 
as measured by DAMA with neutron source integrated over the 6.5 — 97keV and 
the 22 — 330 keV recoil energy range, respectively |1|; Q//) the quenching factors 
evaluated as in Ref. [57] varying as a function of Er] Qm) the quenching factors 
with the same behaviour of Ref. m, but normalized in order to have their mean 
values consistent with Q/ in the energy range considered there. 

Another important effect is the channeling of low energy ions along axes and 
planes of the Nal(Tl) DAMA crystals. This effect can lead to an important devia¬ 
tion, in addition to the other uncertainties discussed above. In fact, the channeling 
effect in crystals implies that a fraction of nuclear recoils are channeled and expe¬ 
rience much larger quenching factors than those derived from neutron calibration 
(see [371 [19] for a discussion of these aspects). Since the channeling effect can¬ 
not be generally pointed out with neutron measurements as already discussed in 
details in Ref. m, only modeling has been produced up to now. In particular, 
the modeling of the channeling effect described by DAMA in Ref. m is able to 
reproduce the recoil spectrum measured at neutron beam by some other groups 
(see Ref. m for details). For completeness, we mention an alternative ehannel- 
ing model, as that of Ref. [88], where larger probabilities of the planar channeling 
are expected. Moreover, we mention the analytic calculation claiming that the 
channeling effect holds for recoils coming from outside a crystal and not from 
recoils produced inside it, due to the blocking effect [89]. Nevertheless, although 
some amount of blocking effect could be present, the precise description of the 
crystal lattice with dopant and trace contaminants is quite difficult and analytical 
calculations require some simplifications which can affect the result. Because of 
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the difficulties of experimental measurements and of theoretical estimate of this 
channeling effect, in the following it will be either included or not in order to give 
idea on the related uncertainty. 

3.4 Migdal effect 

In case of low mass DM particles giving rise to nuclear recoils it is also necessary to 
account for the Migdal effect; this effect is known since long time and is described 
both in dedicated papers [901 E] and in textbooks [92l |93]. It has also been 
recently addressed for the DM field in Ref. [941 [96] . This effect consists in the 

ionization and the excitation of bound atomic electrons induced by the presence 
of a recoiling atomic nucleus. A detailed discussion of its impact in the corollary 
analyses in terms of some DM candidates is given in Ref. |TB]- In this paper, the 
case of mirror nuclei interacting with the target nuclei is considered; thus, since the 
recoiling nucleus can “shake off” some of the atomic electrons, an electromagnetic 
contribution is present together with a recoil signal. Since this contribution is not 
quenched, one can expect that this part (usually unaccounted) can play a role, 
mainly when low mass DM candidates are considered; however, in the present 
case of mirror matter, one can expect second order corrections (of an order not 
exceeding 10“^ on the expected counting rate) when the Migdal effect is accounted 
for. 

3.5 Further uncertainties on parameters 

In the analysis here reported, some discrete cases are considered to account for 
the uncertainties on the measured quenching factors and on the parameters used 
in the nuclear form factors, as already done in previous analyses for other DM 
candidates. The first case (set A) is obtained considering the mean values of the 
parameters of the used nuclear form factors HH and of the quenching factors. 
The set B adopts the same procedure as in Refs. laini, by varying (i) the mean 
values of the measured ^^Na and ^^^I quenching factors up to +2 times the errors; 
(ii) the nuclear radius, and the nuclear surface thickness parameter, s, in the 
form factor from their central values down to -20%. In the last case (set C) the 
Iodine nucleus parameters are fixed at the values of case B, while for the Sodium 
nucleus one considers: (i) ^^Na quenching factor at the lowest value measured in 
literature; (ii) the nuclear radius, rA, and the nuclear surface thickness parameter, 
s, in the SI form factor from their central values up to +20%. 

3.6 Analysis procedures 

As mentioned, the approach exploited by the DAMA experiments is a model- 
independent signature with very peculiar features: the DM annual modulation 
signature. In this case the experimental observable is not - as in other experiments 
- the constant part of the signal, Sq, but its modulation amplitude, Sm, as a 
function of energy. This approach has several advantages; in particular, the only 
background of interest is the one able to mimic the signature, i.e. able to account 
for the whole observed modulation amplitude and to simultaneously satisfy all 
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the many specific peculiarities of this signature. No background of this sort has 
been found or suggested by anyone over more than a decade. Thus, this approach 
does not require any identification of Sq from the total counting rate in order 
to establish the presence of DM particles in the galactic halo. Therefore, the 
DM annual modulation signature allows one to overcome the large uncertainties 
associated to: i) the many data selections/subtractions/statistical-discrimination 
procedures; ii) the modeling of surviving background in keV region; hi) the a priori 
assumption on the nature, interaction type. etc. of the DM particle(s), which are 
necessary in the experiments where the experimental observable is 5o. 

When the DM annual modulation signature is applied, Sq can be worked out - 
for each considered framework ~ by corollary model dependent analysis through a 
maximum likelihood analysis which also takes into account the energy behaviour 
of each detector. However, for simplicity the allowed regions in the parameters 
space (e.g. ^/fe in the DM model described in this paper) can also be calculated 
by comparing - for each k — th energy bin - the measured DM annual modulation 
amplitude with the expectation in the considered framework, In this 

procedure one should remind that the measured counting rate in the cumulative 
energy spectrum - given by the sum of the constant background contribution bk 
and of the constant part of the signal - is about 1 cpd/kg/keV in the lowest 
energy bins; in particular, as discussed e.g. in Ref. m, this constant background 
is estimated to be not lower than ~ 0.75 cpd/kg/keV in the 2-4 keV energy region; 
thus, an upper limit on Sq of ~ 0.25 cpd/kg/keV (5o,max) is derived!^. 

To compare the expectations with the experimental results, the modulation 
amplitudes as function of energy [36] have been considered. The energy bin used 
here is 1 keV and the experimental modulation amplitude in the k — th bin is 
^ We compute the quantity: 


X = 


E 

k 


, ce-xp cth 

y^m.k ‘^1 


m,k) 


+ 


(‘5o,r 


qth \ 
‘^ 0 , 2-41 


® (‘^0^2-4 ~ 'So,max) (13) 


crt 


a. 


2-4 


where the second term encodes the experimental bound about the unmodulated 
part of the signal; here <T 2_4 ~ 10“^ cpd/kg/keV, 0 is the Heaviside function, 
and 5g^2-4 is the average expected signal counting rate in the (2 — 4) keV energy 
interval. The sum in eq. (11311 runs from the software energy threshold (2 keV) to 
20 keV. Given the sharp decreasing shape of the expected signal for the candidate 
considered here, the results are strongly driven by the data points in the (2 — 4) 
keV energy interval. 

The x^ of eq. (|13l) in the Mirror DM model considered here is function of only 
one parameter: v7e; thus, we can define: 

= 0 }. 

The Ax^ is a x^ with one degree of freedom and is used to determine the allowed 
interval of the y/Je parameter at 5 ct from the null signal hypothesis. 

® It is worth noting that not to account for this experimental fact is one of the many reasons 
of incorrect allowed regions put forward by most authors for the particular scenario they adopt. 
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4 Results 


The data have been analysed by taking into account the uncertainties discussed 
in the previous sections; they have been accounted for by evaluating the results 
for the various sets of parameters as summarized in Table [2j In particular, five 


Table 2: Results on the y/Je parameter in the considered scenarios obtained by 
analysing the DAM A data in a mirror DM framework as discussed in the text. 
For each scenario the best fit value of the y/Je parameter and the relative allowed 
interval (corresponding to model providing the deeper Ay^) are reported as well as 
the cumulative allowed interval for y/fe obtained when considering all the above 
mentioned models. The allowed intervals identify the '/Je values corresponding 
to C.L. larger than 5a from the null hypothesis, that is y/Je = 0. See text. 


Scenario 

Quenching 

Factor 

Channeling 

Migdal 

y/Je best 

y/Je interval 
(xlO^®) 

a 

Ql HI 

no 

no 

4.45 X 10^® (9.2(7 C.L.) 

1.86-4.52 
(all) 1.73-114. 

b 

Ql SI 

yes 

no 

2.89 X 10-® (9.3(7 C.L.) 

1.16-2.93 
(all) 0.77-9.72 

c 

Ql SI 

no 

yes 

4.40 X 10-® (9.2(7 C.L.) 

1.85-4.47 
(all) 1.72-107. 

d 

Qii [HZl 

no 

no 

2.44 X 10-® (9.5(7 C.L.) 

1.03-2.48 
(all) 0.94-12.3 

e 

Qlll |S7| -normalized 

no 

no 

5.18 X 10-® (9.0(7 C.L.) 

2.24-5.26 
(all) 1.89-60.1 


scenarios have been considered depending on: i) the adopted quenching factors; 
ii) either inclusion or not of the channeling effect; hi) either inclusion or not of 
the Migdal effect. For each scenario the halo models (138 models as discussed in 
Sect. 13.ip and the relative uncertainties (the three sets described in Sect. 13.Sp 
have been considered. In Table [2] for each scenario the best fit y/fe parameter 
corresponding to the model providing the deeper Ax^ is reported; in addition, 
the allowed intervals of the y/Je parameter for the deeper Ax^ model and for all 
the considered models are reported as well. These allowed intervals identify the 
y/Je values corresponding to C.L. larger than 5a from the null hypothesis, that is 
y/Je = 0. 

In Fig. [2] the distributions of the logiQ{y/Je) allowed intervals of all the models 
are shown for each considered scenario. The scenarios a, c and e are very similar, 
while the scenarios either with channeling effect (b) or with with the Qn quenching 
factors (d) support lower values of the Mirror DM coupling. 

In Fig. [3] comparisons between the DAM A experimental modulation ampli¬ 
tudes and some expectations for Mirror DM are shown. 

It is worth noting that in all the considered scenarios for mirror DM the DAMA 
signal in the 2-6 keV energy interval arises from interactions mainly with Sodium 
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Figure 2: Distributions of the logiQ(v7£) intervals obtained for all the models in 
each considered scenario, identified by the letter, according to Table [2j See text. 


nuclei 0. This effect is due to the fact that the considered Mirror DM particle is 
quite light: Ma' — 5mp. 

The cumulative allowed intervals of the y/Je parameter selected by the DAMA 
data for each scenario (see Table [2]) are depicted in Fig. [H where also the over¬ 
all allowed band is shown. The obtained values of the y/Je parameter are well 
compatible with cosmological bounds cited in the introduction. 

Finally, we would like to note that the highest C.L. among all the analysed 
cases (see Table ED is obtained for the scenario d where: i) the quenching factors 
values are according to Ref. m; ii) the halo model is the no-rotating Evans’ 
logarithmic C2 model (see Tabled]) with vq = 270 km/s, pdm = = 1-68 GeV 

cm“^ and Vesc = 650 km/s; iii) set B, as defined in Sect. 13.51 

If the assumption Ma' — 5mp is released, the allowed regions for the y/Je 
parameter as function of Ma' can be obtained by marginalizing all the models for 
each considered scenario as given in Table[2j This is shown in Fig. [5|where the Ma' 
interval from few GeV up to 50 GeV is explored. These allowed intervals identify 
the y/Je values corresponding to G.L. larger than 5a from the null hypothesis, that 
is y/Je = 0. The five scenarios defined in Tabled) can be recognized on the basis 
of different hatching of the allowed regions; the black line is the overall boundary. 

^For example, the lodine/Sodium contribution ratio in the best fit case of the scenario d is 
0.0062 and 2.4 x 10“® for the first two energy bins, (2-3) and (3-4) keV, respectively. 
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Figure 3: Examples of expected modulation amplitude for the Mirror DM can¬ 
didate considered here. In particular, there are shown the best fit cases of the 
five scenarios (from top to bottom): i) Mirror DM candidate with Qi quenching 
factors and without channeling effect (scenario a); ii) Mirror DM candidate with 
Qj quenching factors and channeling effect included (scenario 6); iii) Mirror DM 
candidate with Migdal effect included in the interaction (scenario c); iv) Mirror 
DM candidate with Qn quenching factors (scenario d); v) Mirror DM candidate 
with Qiii quenching factors (scenario e). 
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Figure 4: DAMA allowed intervals for the y/Je parameter, obtained by marginal¬ 
izing all the models for each considered scenario as given in Table [2j The overall 
range is also reported. 

5 Conclusions 

The model-independent annual modulation effect measured by the DAMA Collab¬ 
oration, which fulfills all the requirements of the DM annual modulation signature, 
has been examined in the context of asymmetric mirror matter model interact¬ 
ing with the ordinary nuclei via the photon-mirror photon kinetic mixing portal, 
1 We have assumed that mirror atoms constitute a fraction / of the DM 

in the Galaxy, and we have derived the allowed physical intervals for the combina¬ 
tion of parameters y/fe, accounting also for various of the existing uncertainties. 

The allowed values for y/ye in the case of mirror hydrogen atom, Z' = 1, 
ranges between 7.7 x to 1.1 x 10“^. The values within this overall range 

are well compatible with cosmological bounds. In particular, the best fit values 
among all the considered scenarios gives vT^b.f. = 2.4 x 10“®. If the assumption 
Mj^i ~ 5 GeV is relaxed, the allowed regions for the y/Je parameter as function 
of Mjyi have been obtained by marginalizing all the models for each considered 
scenario. We have also to remark that the atomic form-factor of mirror was not 
taken into account which is correct if mirror electron is light enough. In the case it 
is heavy, the obtained results for y/Je should be rescaled by the square root of the 
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Figure 5: Allowed regions for the ^/Je parameter as function of when the 
assumption M^' — 5mp is released, obtained by marginalizing all the models 
for each considered scenario as given in Table [5J The M^i interval from few 
GeV up to 50 GeV is explored. These allowed intervals identify the ^/Je values 
corresponding to G.L. larger than 5a from the null hypothesis, that is y/fe = 0. 
The hve scenarios dehned in Table [2] can be recognized on the basis of different 
hatching of the allowed regions; the black line is the overall boundary. 


mirror atomic form-factor. In the case of helium-like A-atoms with Z' = 2, the 
ranges of y/Je obtained for the hydrogen case should be rescaled up by a factor 2. 

In our consideration, we did not take into account the possible modifications in 
the dark matter distribution in the Galaxy due to self-scattering properties of dark 
matter. However, marginalizing over a vast amount of models that we considered 
in our conservative approach should supersede all possible uncertainties that can 
be induced by this reason. Let us remark also that the photon-mirror photon 
kinetic mixing, besides the possibility of dark matter direct detection, could bring 
to other interesting cosmological consequences. In particular, as it was shown in 
Ref. [98], in rotating protogalaxies the circular currents could emerge due the 
Rutherford-like scattering of the electrons with dark mirror particles in the dark 
matter halos, with cross-sections suppressed by a parameter e. After a moderate 
dynamo amplihcation, these currents could give origin to the observed galactic 
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magnetic fields of few /xG with the coherence scales of order of 1 kpc. Photon- 
mirror photon kinetic mixing portal and search for mini-charged particles in whole 
is elusive for LHC and other precision experiments are needed for detecting their 
production in laboratory conditions. 
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